Confined chaotic behavior in collective motion for populations of 

globally coupled chaotic elements 



Naoko NAKAGAWA 

Laboratory for Information Synthesis, RIKEN Brain Science Institute, Wako, Saitama 351-01, 

Japan 

OO 

Teruhisa S. KOMATSU 

^ , 

^ , Department of Pure and Applied Sciences, College of Arts and Sciences, University of Tokyo, 

, Tokyo 153, Japan 



O ■ Abstract 

\o 

O 

00 . 

, The Lyapunov exponent for collective motion is defined in order to charac- 

. 

terize chaotic properties of collective motion for large populations of chaotic 
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collective motion is found to be confined within a small scale, whose size is 

estimated using the value of the Lyapunov exponent. Finally, we conjecture 

why the collective motion appears low dimensional despite the actual high 

dimensionality of the dynamics. 
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I. INTRODUCTION 



Systems consisting of large populations of interacting dynamical elements are widely 
distributed in nature, from communities of ants to biological cell assemblies and neural 
networks. Many such populations are similar in the respect that they exhibit various kinds 
of collective behavior. One possible approach to the mathematical study of such collective 
behavior is to concentrate on certain idealized models such as interacting limit cycles and 
chaotic maps. An abundance of literature has been devoted to the study of these models 
|T|-0. It is suggested in some foregoing studies |]^-^ that a population of this type as a 
whole exhibits low dimensional behavior. This seems to be true even when the individual 
elements appears to be mutually uncorrelated, a situation in which the population could only 
be regarded as a dynamical system of extremely high dimension. It still remains unclear 
whether such collective behavior could be understood in terms of low-dimensional dynamical 
systems. 

In globally coupled tent maps, collective motion has been understood partly from a 
macroscopic viewpoint In these studies, the collective motion was discussed only 

with regard to its size, without consideration of its detailed structure. Some beautiful 
relations regarding the macroscopic properties of the collective motion were explored and 
a phase diagram was proposed. These results, however, do not represent progress toward 
an understanding of collective motion in terms of low- dimensional dynamical systems. In 
this Paper, we first focus on the dynamical properties of collective motion. We define the 
Lyapunov exponent characterizing this collective motion and then a special relation involving 
this quantity is found. This relation suggests that the collective motion is always chaotic, 
whenever it appears. Furthermore, chaotic behavior in the collective motion are found to be 
confined within a small scale whose size is related to the value of the Lyapunov exponent. 
Finally, on the basis of our numerical results, we conjecture an answer to the question: Why 
does the collective motion appear low dimensional despite the actual high dimensionality of 
the dynamics? 
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II. MODEL : GLOBALLY COUPLED TENT MAPS 



As model systems, we have chosen globally coupled tent maps, because these systems 
are particularly well suited for detailed numerical analysis: We can deal with an ideal limit 
of infinitely large population. 

Globally coupled maps (GCM) are given by an assembly of elements whose behavior 
is determined by N identical maps with all-to-all coupling. The individual elements are 
then under the influence of a common internal field which may be referred to as a mean field 
0. We assume that a single isolated element evolves according to Xn+i = f{Xn), where 
n designates discrete time steps. Under the interaction through the mean field the i-th 
element is then assumed to evolve as: 

xi:i = fix^:^)+KK, (1) 

where K is the coupling strength. In this Paper, we consider the situation in which / is a 
tent map: 

f{X) = -a\X\ + ^, (2) 
and the mean field hn is defined as: 

1 ^ 

hn^-Y^fiXj^-^). (3) 
i=i 

Thus, our system is characterized by the two parameters a and K in addition to the total 
number of elements A^. Each tent map has a band splitting point a = \/2 and we thus 
stipulate that a satisfies < a < 2. It is thereby ensured that the population will never 
split into sub-populations. If the system size N is finite, this finiteness becomes the source of 
fluctuations of the mean field. Such an effect may obscure pure collective motion. Thus we 
work with the limit of large A^. In fact, we confirmed that finite size effects can be regarded 
as noise acting on the pure collective motion. For — oo, the population dynamics of 



GCM can be described by the Frobenius-Perron equation ||TT| for the distribution p{X): 
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Pn+l{X) = JSiX- fix') - KK)pn{X')dX' 

K = j fix')prXx')dx'. 



(4) 



(5) 



We worked out a numerical scheme for the exact integration of Eqs.(^ and (|^) whose 
precision is hmited only by round-off errors [|T^. Numerical results in this Paper were 
obtained using only a single initial distribution, uniform over the interval [/ o /(O) : /(O)]. 

III. GROWTH OF COMPLEXITY IN COLLECTIVE MOTION 

Collective motion can be observed by studying the dynamics of the order parameter hn 
given by Eq.(§). Figures ^a) and |l](b) are return maps for this order parameter. There 
we see that the fluctuations of /i„ undergo quasi-periodic motion for small K, but for larger 
K they display more complicated motion. Roughly speaking, larger values of K generically 
result in more complicated collective behavior over all values of a. 

In the case of small K, the collective motion appears quasi-periodic. Despite the familiar 
appearance of this collective motion, we should discuss carefully the nature of these dynam- 
ics, because the system has an infinitely high dimension. The chaotic structure on the torus 
seems to grow gradually, without an indication of bifurcation from quasi-periodic to the 
chaotic behavior. It may be the case that the structure on the torus persists as long as there 
is a non-vanishing interaction. In following sections, we focus on the chaotic properties of 
the collective motion. 



In order to characterize the dynamical properties of the collective motion, we wish to 
define a value similar to the Lyapunov exponent to characterize the collective motion. We 
note that it is unfeasible to calculate conventional Lyapunov exponents based on the phase 
space structure of the system, because our system is of an infinitely high dimension. However, 
what we wish to concentrate on is not the phase space structure but the dynamical properties 



IV. EXPONENTIAL GROWTH OF DISTURBANCES 
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of the collective motion. Hence we attempt to define a quantity similar to the Lyapunov 
exponent to describe the dynamics of the collective motion. 

We first add a small disturbance to the distribution Pn{X) at time step n. The mean field 
hn is changed to /i^'^-* as a result. The difference between the disturbed and non-disturbed 
mean fields, 

di = \hl^li-hr,+i\, (6) 

depends on I, the time lapsed since the application of the disturbance. Without an inter- 
action {K = 0) among the elements in the population, the distribution p''^]^i{X) gradually 
approaches that of the non-disturbed case pn+i{X). This is because the two distributions 
should realize the same invariant measure of a single tent map. Thus — * as / becomes 
large. When K ^ 0, collective motion generally occurs. We expect that its dynamical 
properties can be represented by the behavior of 6i. If the collective motion is completely 
quasi-periodic, the difference 6i will remain on the same order as 6q, while, if it is chaotic, 
6i will grow exponentially with the increase of /. 

In numerical experiments, 6i is estimated for one set of parameters (a, K) as follows: 
The value \h^n+i ~ ^n+i\ is calculated from various initial distributions of p';^\X) with a fixed 
small difference 6o = \h!;f^ — hn\. The sample average over various \h^n+i ~ hn+i\ calculated 
in this way gives 6i. In our computations, this average was calculated from either 10 or 
100 initial disturbed distributions. We point out here that individual l^^'!^; — hn+i\ actually 
depends on / in a manner similar to 6i. 

In Fig.^ we display examples for the exponential growth of the difference 6i between the 
mean fields. The initial difference 6o is fixed at as small as 10^^°. The difference Si exhibits 
a clear dependence on the time step /: 6i grows exponentially and finally saturates near a 
certain scale. It is surprising that the exponential growth of 6i is observed for a considerably 
large interval of /. It is therefore possible to define the exponent A which characterizes the 
exponential growth of 6i. The exponent A corresponds to the initial slope in Fig.|^. For 
values of / at which 6i displays exponential growth, we can define A according to 
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A^Mf). (7) 

This exponent was found to be independent of the initial disturbance of the distribution 
for sufficiently small values of 6q and sufficiently large values of /. Hereafter, we call the 
exponent A the Lyapunov exponent for the collective motion. 

In Fig.0 we display three examples of 6i, where the values of K are 0.6, 0.8 and 1.0 while 
the value of a is fixed according to — log(j(2 — a) = 4. Larger values of K results in larger 
values of A. The collective motion is concluded to be of a chaotic nature in every case of 
Fig.^]. If we were to calculate conventional Lyapunov exponents based on the phase space 
structure of this system, a maximal exponent with a value close to log a, the Lyapunov 
exponent for a single tent map, would be obtained. We note that the value of A is generally 
smaller than log a. In the case of Fig.|^, log a = 0.48, which is larger than the values found 
for A (for instance A = 0.17 when K = 1.0). This implies that the dynamics of the collective 
motion are independent of the dynamics of the individual elements. 



V. HILLY STRUCTURE OF THE LYAPUNOV EXPONENT 

Figure ^ shows the Lyapunov exponent A for the collective motion as functions of a ( or, 
more precisely, as functions of m = — log^(2 — a) ) at fixed values of K. Here, hilly structure 
is displayed by the Lyapunov exponent A, as seen in Fig.^(a). For comparison, we display the 
hilly structure corresponding to macroscopic properties of the collective motion in Fig.^(b). 
These macroscopic properties are represented by the amplitude F of the collective motion, 

F ^ ^{{K-{K)Y). (8) 

Here () represents a long time average. It is worthy to note that the peaks and valleys of 
these two sets of hilly structure seem to be located at the same values of a. For further 
detail concerning the hilly structure of F, see |p. Two special series of parameter values of 
a, the 'golden' and 'silver' values, are the key to understanding the hilly structure. Briefiy, 
the golden values of the parameter a are situated in the middle of each hill and the silver 
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values are at the minimum points in the valleys. The same situation is found in the case of 
the Lyapunov exponent A for the collective motion. 

We now concentrate on the differences between these two sets of hilly structure. The 
hilly structure of the Lyapunov exponent A exhibits gentle ups and downs over the range of 
values of a, while the amplitude F displays very sharp ups and downs. The most important 
difference appears in the valleys of the hilly structure: The value of F decays linearly to 
in each valley, but the value of A in each valley seems to remain on the same order as that 
assumed at the top of the nearest hill. However, we will later show that the values of A 
in fact also decay to 0. Decaying will be found to occur only in a very narrow region of a 
around the minimum point in each valley. 

The regions in which the value of A is close to zero are sufficiently narrow that we can 
ignore them in estimating the order of the Lyapunov exponent A as a function of a. For 
fixed values of K the value of the Lyapunov exponent A decays as a becomes large. With 
the definition m = — log„(2 — a), this decaying tendency is estimated as 

A ~ (9) 

for positive values of K. In Section VII, we discuss the exceptional regions of a (i.e., those 
in which this relationship does not hold). 

VI. SCALING RELATION OF THE LYAPUNOV EXPONENT 

In this section, we discuss the relation of the exponent A with the coupling strength K. 
The value of a is confined to be one of the golden values, while the value of K is changed. 
Here a golden value, say gp, is a value of a such that one isolated tent maps possesses the 
property that a trajectory beginning at the peak C of the tent map returns to C after p 
steps. There exists a sequence of such golden values, parameterized by p. As shown in 
Fig.^, it is confirmed that the collective motion becomes increasingly chaotic as the value 
of K increases. For all the golden values we considered, we have found that the Lyapunov 
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exponent A grows with the couphng strength K according to the relation A ~ . The 
value of A also obeys the relation Eq.(H). Thus we find that the Lyapunov exponent A for 
the collective motion behaves as 

A ~ — , (10) 

m 

for the golden values of a. While we have found this relation by confining our consideration 
to these golden values, it should be noted that a relation similar to Eq.(|lOD is expected to 
hold for all values of a except in very narrow regions, due to the gentleness of the hilly 
structure of A. Also note that Eq.(0) in some sense describes the envelope of A in the 
(a, K) parameter space, as it does not include the periodic dependence on m or the fine 
peak structure of A. 

Equation ( [T0| ) indicates that when there is an interaction between elements (i.e. X 7^ 0) 
the collective motion is always chaotic. Thus collective motion displaying quasi-periodicity, 
as is seen in Fig.|l|(a), is concluded to actually be chaotic. As the coupling strength K 
approaches or the parameter a approaches 2 (i.e. m — > 00), the Lyapunov exponent A for 
the collective motion vanishes asymptotically. 

VII. LYAPUNOV EXPONENTS AROUND DISAPPEARING POINTS OF THE 

COLLECTIVE MOTION 

The scaling relation for the Lyapunov exponent A, Eq . (Il0|) , leads us to conclude that 
the collective motion is always chaotic. However, there do exist exceptional cases in which 
the collective motion disappears. In particular, it disappears at the minimum points in 
each valley of the hilly structure of the amplitude F. Around these minimum points, the 
Lyapunov exponent for the collective motion is expected to decay to vanishing or negative 
values. However, it is difficult to detect this behavior in Fig.^. We now investigate in detail 
the nature of the Lyapunov exponent around the minimum points. 

The values of a at the minimum points are called 'silver values'. Here, silver values 
constitute another particular series of parameter values of a. Figure ^ shows the decaying 
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behavior of the Lyapunov exponent around a silver value s„ of the parameter a. \Aa 



neighborhood of s„. The value of A becomes small as |Aa| approaches 0. As reflected by 
this figure, the Lyapunov exponent scales with the distance |Aa| as 



This relation shows that the Lyapunov exponent decays to as the value of a approach to 
a silver value. The Lyapunov exponent is non-positive only at the silver value. There, the 
collective motion does not appear, so that chaotic nature of this motion is never realized. 
What we wish to emphasize is that A falls to zero in only a very narrow neighborhood of the 
silver value (see Fig.^(b)). For all other values of a, the Lyapunov exponent is approximately 
given by Eq. (|1^) . The narrowness of the decaying region around each silver value causes the 
apparent smoothness of the Lyapunov exponent as a function of a. The motion of the system 
is chaotic everywhere, except at the minimum point of each valley, where collective motion 
does not exist. We now conclude that, whenever collective motion appears (i.e.F 7^ 0), it is 
chaotic. 

VIII. MINIMUM STRUCTURE OF THE COLLECTIVE MOTION 

Let us now turn in our investigation of the Lyapunov exponent for the collective motion 
to consideration of the characteristic scale, that is, the scale dc at which the exponential 
growth of 61 saturates (see Fig.|^). The value of dc is determined numerically through the 
time average of Si over an interval of sufficient length just after the saturation. After the 
saturation of 61, it may continue to increase very slowly, until it becomes comparable to 
the size F of the collective motion. We, however, can choose an appropriate interval for 
obtaining a valid value of dc through the time average of 61 just after saturation. 

Below the scale characterized by 6c, the collective motion behaves chaotically correspond- 
ing to a linear instability of the collective motion. The structure characterized by 6c is of a 



represents the distance from a to the silver value 



Sn- Figure ^ displays a very narrow 



A ~ (—log |Aa|) 3. 



(11) 
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scale at which the nonhnearity first appears in the collective motion. The chaotic behavior 
of the collective motion is confined in this structure. We hereafter call the structure of the 
collective motion characterized by 6c the minimum structure. We determined an expression 
for the characteristic scale dc of the minimum structure. 6c was found to depend on the 
value of K according to a certain scaling relation. As is seen in Fig.^ for sufficiently small 



values of K, 6c behaves as 



4~exp(--i^). (12) 



This figure was obtained at a particular golden value of a while the relation (|l^) was con- 
firmed at every golden value we considered. The relation in Eq.(0) indicates that the 
minimum structure exists however small the coupling strength K is. 

Let us now compare the size 6c of the minimum structure with that of the collective 
motion. The amplitude F of the collective motion explored in obeys a different type of 
scaling relation, 

i^~exp(-l). (13) 

The scaling relations for both 6c and F are given by monotonically increasing functions of 
K. Both sizes decay to as decreases. However, the rate at which this decay occurs is 
much different for 6c and F. The size 6c of the minimum structure becomes much smaller 
than the amplitude F of the collective motion for sufficiently small values of K, while they 
are comparable for if ^ 1. The minimum structure and the chaotic behavior within it 
belong to a collective set of phenomena in the submacroscopic scale: They are properties 
characterizing the smaller scale of the collective motion, while the collective motion is a 
product of the macroscopic nature of the distribution p„(X). 

IX. LYAPUNOV EXPONENT AND THE MINIMUM STRUCTURE 



For the golden values of the parameter a the size of the minimum structure is found to 
be strongly correlated with the value of the Lyapunov exponent A for the collective motion: 
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5e ~ exp(--). (14) 

This relation is obtained through the substitution of Eq.(|IO|) into Eq.(]T^). We conjecture 
that this relation is not mere coincidence (i.e., its validity is not confined to the golden 
values) but holds for all sets of (a, K). The characteristic size 6c corresponding to the chaotic 
behavior should naturally be related to the Lyapunov exponent. This conjecture has been 
confirmed through numerical comparison of A and — (log^ 4)~^- In Fig.^ we show that these 
two values depend on a in a very similar manner, although the values of — (log^^c)"^ are 
more scattered than those of A. We thus conclude that the relation 

6c ^ (15) 

is valid for any set of (a, K). 

The relation in Eg. ([T5|) implies some properties of the minimum structure. As was ar- 
gued, the Lyapunov exponent A remains positive except at the silver values of a. Accordingly, 
the minimum structure exists for almost all sets of {a,K), where it becomes small as the 
chaotic properties become weak. When the coupling strength K approaches 0, the values of 
the Lyapunov exponent was found to approach with the relation in Eq.(|l^). Correspond- 
ing to this, the size 6c of the minimum structure also decreases to with a much larger 
decreasing rate with K, as was described in Eg. ([T^) . At the silver values of a, the minimum 
structure does not exist {6c = 0), because the collective motion disappears at these points. 
The value of A at these points was numerically found to be 0, which is consistent with the 
non-existence of the minimum structure {6c = 0). The relation in Eq.(|15D indicates further- 
more that the size 6c of the minimum structure also exhibits hilly structure as a function of 
a, like A and F. In fact, the hilly structure of 6c has been numerically generated. However, 
because the values of 6c found in this study are widely scattered, the functional dependence 
of 6c on a was not made completely clear. 

It was concluded that the collective motion is almost always chaotic even though it 
appears quasi-periodic. We now obtain a second conclusion from the relation in Eq.(|l5D 
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that the chaotic behavior is always confined in a submacroscopic scale given by Eq.(^). 
If the coupling strength K is sufficiently small, the chaotic behavior is observed only in a 
much smaller scale than the scale of the collective motion. We argue in a following section 
that this confinement of the chaotic behavior leads to the low dimensional appearance of 
the collective motion. 

X. CHAOTIC FLUCTUATIONS IN THE COLLECTIVE MOTION 

The minimum structure is related to a certain observed structure in the collective motion. 
Consider first the case that the collective motion locks to some periodic states, as is seen 
in Fig.^]. Each periodic point possesses a region which is subject to strong fiuctuations as 
is shown in the blowups. The extent of these fiuctuations seems to be comparable to the 
scale dc of the minimum structure. When, as in Fig.|I], the collective motion does not lock to 
periodic states, the structure corresponding to 6c is similar to the case of the locked periodic 
states. As is typically observed (see Fig.|l|(b)), the collective motion fiuctuates around some 
quasi-periodic motion. The scale of these chaotic fiuctuations seems to be on the same order 
as the value of Sc- 

Accepting the correspondence between the extent of the chaotic fiuctuations and the 
size of the minimum structure, we can interpret the observed collective motion in terms of 
the macroscopic properties of the collective motion and the minimum structure. If K > 1, 
dc is comparable to the size F of the collective motion. In this case the collective motion 
becomes fully-developed chaotic motion (see Fig.|l](b)). When K <^ 1, the size dc of the 
minimum structure becomes much smaller than the size F of the collective motion. Here, it 
is possible for the collective motion to approach a form similar to that displayed by a low- 
dimensional system, as in the case of the quasi-periodic motion depicted in Fig.|l|(a). Our 
results imply that it is not possible for the system to realize quasi-periodic behavior which 
lacks chaotic fiuctuations. Whenever the collective motion appears, it is always accompanied 
by chaotic fiuctuations. It is important that these chaotic fiuctuations are confined within 
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the submacroscopic scale whenever they appear. The collective motion does not realize true 
quasi-periodic behavior, however, it reduces to a form which appears to be quasi-periodic due 
to the confinement of the chaotic fiuctuations within the submacroscopic scale. As a result, 
the quasi-periodic behavior of the collective motion becomes conspicuous, especially for 
small values of K, where the characteristic scales of the collective motion and the minimum 
structure are very different. 

We now make a conjecture to answer the question of why the collective motion appears 
low dimensional despite the actual high dimensionality of the dynamics. Recall now that 
the collective motion results from dynamics of infinitely high dimension. This is because the 
present model (Eq.(^ and (|^)) consists of an infinite number of interacting chaotic elements 
without any mechanism to reduce the number of the degrees of freedom. It is natural to 
imagine that the chaotic fiuctuations in the submacroscopic scale constitute chaotic behavior 
of infinitely high dimension. Due to the confinement of the high dimensional chaos within 
the submacroscopic scale, the collective motion reduces to a form which appears to be low- 
dimensional. 

XI. MACROSCOPIC BEHAVIOR VS CHAOTIC FLUCTUATIONS 

We now use a term macroscopic behavior definitely for the coarse-grained behavior of 
the collective motion below the submacroscopic scale, i.e. below the scale of 6c of the 
minimum structure. Macroscopic behavior was observed to exhibit low-dimensional motion 
such as quasi-periodic and locked periodic motion. We now briefiy discuss the variety of the 
macroscopic behavior for the collective motion. The macroscopic behavior of the collective 
motion changes very frequently in the parameter space {a,K): typically the quasi-periodic 
motion repeatedly locks to and releases from some periodic states (see Fig.^). Examples of 
the collective motion are displayed in Fig.^. The locking phenomena of the collective motion 
are often followed by the instability of the locked points. We find an assembly of tori, as in 
Fig.^(a). These tori can also become locked, where each torus in the assembly turns into 
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to locked states so that an assembly of the locked periodic states is formed. In such a way, 
a kind of hierarchical structure of locked states can be realized when the creation of tori, 
resulting from the destabilization of periodic points, is repeated again and again, leading to 
the formation of assemblies of tori. When the collective motion does not lock, it comes to 
possess fine structure reflecting a delicate dynamical property for the macroscopic behavior 
(Fig.|(b)). 

As is seen in the above examples, the macroscopic motion can become complicated, 
in some cases possessing fine structure, such as a kind of hierarchical structure. We here 
emphasize that, without regard to the dynamical properties of the macroscopic scale, the 
minimum structure inevitably exists along with the macroscopic behavior of the collective 
motion. For all patterns of the macroscopic behavior, the chaotic behavior only remains 
below the submacroscopic scale. In this scale the dynamical properties of the macroscopic 
behavior is destroyed, whereas the chaotic properties in the submacroscopic scale do not 
seem to affect the collective behavior in the macroscopic scale. The appearance of the 
minimum structure marks the disappearance of the detailed structure of the macroscopic 
behavior. The existence of the submacroscopic scale reminds us of the Kolmogorov scale 
in turbulence, where the structure of the energy cascade is destroyed by the mechanism of 
thermal dissipation. 

XII. SUMMARY AND DISCUSSION 

In this paper we discussed the collective motion of globally coupled tent maps from the 
viewpoint of the dynamical properties of the distribution. We defined the Lyapunov expo- 
nent for the collective motion from the growth of a disturbance in the distribution Pn{X) 
through projection on the order parameter /i„. The values obtained for the Lyapunov expo- 
nent are nontrivial in the sense that they cannot be inferred from the dynamical properties 
of one tent map. It was found that the collective motion is chaotic whenever it appears, 
even though it seems to be periodic or quasi-periodic. The chaotic behavior of the coUec- 
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tive motion seems to exist inside the minimum structure of the collective motion, where 
disturbances grow exponentially. This minimum structure persists for any collective motion 
with chaotic behavior and its size can be estimated from the value of the Lyapunov expo- 
nent. These results suggest the reason that the collective motion exhibits low-dimensional 
behavior despite the actual high dimensionality. 

It is necessary to explain how the minimum structure come into existence. The scaling 
relations for the size 5c of the minimum structure, Eg . ([12D , and the Lyapunov exponent A 
for the collective motion, Eq. ([lOD , are similar to the forms obtained analytically by Ershov 
and Potapov ||13|. They considered the size of the fluctuations and Lyapunov exponents for 
globally coupled tent maps. Due to their abstruse analysis, we are not sure if the fluctuations 
considered there are related to the minimum structure discussed in this Paper. We need 
a more intuitive explanation regarding the occurrence of the chaotic fluctuations and their 
confinement within the minimum structure. 

This research is partially supported by JSPS Research Fellowships. 
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FIGURES 

FIG. 1. Two types of collective motion at a = 55 (defined by — log(j(2 — a) = 4) for the 
return map of hn- (a) For K = 0.1, we find a torus representing quasi-periodic motion, (b) For 
K = 1.0, the return map displays more complicated torus-like motion, possibly accompanied by 
fine structure. 

FIG. 2. Exponential growth of the difference 61 between the order parameters h^^i and 
The vertical scale is logarithmic. Here a = (— \ogg^{2 — a) = 4), while K = 0.6, 0.8 and 1.0. All 
points represent sample averages of 10 different initial distributions {X) with a fixed difference 
So = 10-30. 

FIG. 3. Hilly structure as a function of m for if = 1.0 and 0.5 (x). Here m = — log„(2 — a), 
(a) m vs Lyapunov exponent A. In the insert, the same data are displayed on logarithmic scales 
for both axes. The dotted line indicates A = m~^. (b) m vs amplitude F of collective motion. 
There are several hills and valleys whose locations seem to correspond to those in (a). Some of the 
golden values of a coincide with integer values of m while some of the silver values coincide with 
the minimum points of valleys of F. 

FIG. 4. Scaling relation of the Lyapunov exponent A for the collective motion with K. Here A 
vs K is shown for three golden values gp, p = 3 - 5. a = gp is equal to — log„(2 — a) = p — 1. The 
dotted line indicates A = K^. 

FIG. 5. Decay of the Lyapunov exponent A for the collective motion near a silver value 
sj = 1.81327 • • •. Here Aa = a — si. (a) — log | Aa| vs A with logarithmic scales for both axes. The 
dotted line corresponds to A oc (— /o5r|Aa|)-3 . (b) |Aa| vs A with normal scales for both axes. The 
region where A decays to is found to be very narrow when A is plotted in this way. 

FIG. 6. Scaling relation of the size 6^ of the minimum structure with K. Here 5c vs 1/if^ is 

shown for the golden value a = g^ (—\ogg^{2 — a) = 3). The vertical scale is logarithmic. The 
dotted line indicates dc oc exp{—K~'^). 
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FIG. 7. Relation between the Lyapunov exponent A and the size 5c of the minimum structure. 
Here m = — log„(2 — a). This figure is for the case K = 0.5. (a) The Lyapunov exponent A for the 
cohective motion, (b) A' = (— 2 log^ (5c)~^, where 6c is the size of the minimum structure (also see 
Figj2|). We find good correspondence between (a) and (b). 

FIG. 8. Locking phenomenon of the collective motion at a = (74 (\og^{2 — a) = 3) and K = 0.5, 
for the return map of In the blowups, fluctuations are found around each locked point. The 
extent of the fluctuations is comparable to the size 6c of the minimum structure at a corresponding 
set of parameters {a,K). 

FIG. 9. Examples of collective motion. Here a = §4 (— log(j(2 — a) = 3). Each figure is for 
the return map of (a) K = 0.501. Each periodic point in Fig.^ becomes a torus due to the 
instability of the locked point. For a slightly larger value of K, these small tori exhibit period 
locking again, (b) K = 0.507. The locked state is not stable here. A kind of delicate structure 
is observed, reflecting the stability of the macroscopic dynamics. Points are seen to be scattered 
when a pattern is observed at a sufficiently small scale. These fluctuations appear even after a 
transient behavior has disappeared. 
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